Phase diagram of the random frequency oscillator: 
The case of Ornstein-Uhlenbeck noise 
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We study the stability of a stochastic oscillator whose frequency is a random process with finite 
time memory represented by an Ornstein-Uhlenbeck noise. This system undergoes a noise-induced 
bifurcation when the amplitude of the noise grows. The critical curve, that separates the absorbing 
. . . phase from an extended non-equilibrium steady state, corresponds to the vanishing of the Lyapunov 

exponent that measures the asymptotic logarithmic growth rate of the energy. We derive various 
\ expressions for this Lyapunov exponent by using different approximation schemes. This allows us 

. to study quantitatively the phase diagram of the random parametric oscillator. 
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I \\ I. INTRODUCTION 

of the control parameter, the critical value of the threshold can change and noise can delay or favor a phase transition 
0. In the first case, randomness can be useful as a stabilizing mean and in the second, noise may help trigger a phase 
I transition that is otherwise very difficult to achieve; for example, in the dynamo effect, the role of the noise generated 
by fluid turbulence is not well understood at present and it is possible that the critical magnetic Reynolds number 
decreases with noise 0- In certain cases, a physical system subject to noise undergoes bifurcations into states that 
, have no deterministic counterparts: the stochastic phases generated by randomness have specific characteristics (such 
as scaling behavior or critical exponents) that define new universality classes Q ■ 
r \ One of the simplest systems that can be used as a paradigm for the study of noise-induced phase transitions is 
■ O the random frequency oscillator For instance, in practical engineering problems, the Duffing oscillator with 

^ \ random frequency has been used as a model to study stability of structures subject to random external forces, such as 
2 . earthquakes, wind or ocean waves |lOl 111. 12, 13] . Whereas a deterministic oscillator with damping evolves towards 
I— — i' the unique equilibrium state of minimal energy, the behavior changes if the frequency of the oscillator is a time- 
dependent variable. Due to continuous energy injection into the system through the frequency variations, the system 
may sustain non-zero oscillations even in the long time limit. The case when the frequency is a periodic function of 
time is the classical problem of parametric resonance known as the Mathieu oscillator; the phase diagram is obtained 

■ by calculating the Floquet exponents defined as the characteristic growth rates of the amplitude of the system [l3 |. 
\ When the frequency of t he p endulum is a random process, the role of the Floquet exponents is taken over by the 

00 Lyapunov exponents flS^, 'l6'|. The system undergoes a bifurcation when the largest Lyapunov exponent, defined as 
the growth rate of the logarithm of the energy, changes its sign. Thus, the Lyapunov exponent vanishes on the critical 
^ : surface that separates the phases in the parameter space. This criterion involving the sign of the Lyapunov exponent 
has a firm mathematical basis and clarifies the ambiguities that were found in the study of the stability of higher 
moments [isllr^ . 

In a recent work [l^ , we have carried out an analytic study of the phase diagram of the random oscillator driven 
by a Gaussian white noise frequency. We have shown that, in the case of an inverted pendulum, the unstable 
• fixed point can be stabilized by noise and a noise-induced reentrant transition occurs. These results are based on an 

■ exact formula for the Lyapunov exponent poL l2ll |23 | . In the present work, we intend to study the phase diagram 
O I of an oscillator whose frequency is a random process with finite time memory. More precisely, we consider here the 
J-^ ' case of an Ornstein-Uhlenbeck noise of correlation time r. From a physical point of view, the influence of a finite 
^ , correlation time on the phase diagram is an interesting open question : does a finite correlation time favor or hinder a 

noise-induced transition? In particular, we wish to determine how the shape of the transition curve is modified when 
the noise is colored. In the white noise case, the asymptotic behavior of the critical curve when the amplitude of 
noise is either very small or very large is known explicitly and presents a simple scaling behavior 18]. How do these 
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scalings change when the noise is correlated in time? 

Due to the finite correlation time of the noise, the random oscillator is a non-Markovian random process and there 
exists no closed Fokker-Planck equation that describes the dynamics of the Probability Distribution Function (P.D.F.) 
in the phase space. This non-Markovian feature hinders an exact solution in contrast with the white noise case where 
a closed formula for the Lyapunov exponent was found. We shall therefore have to rely on various approximations 
to carry out an analytical study of the phase diagram. The results obtained by different approximations will be 
compared with numerical results and with an exact small noise perturbative expansion. The various approximations 
have different regions of validity in the parameter space : this allows us to derive a fairly complete picture of the 
phase diagram of the random oscillator subject to an Ornstein-Uhlenbeck multiplicative noise. 

The outline of this work is as follows. In section 2, we derive general results about the stochastic oscillator with 
random frequency: thanks to dimensional analysis, we reduce the dimension of the parameter space from four to two 
and show how the Lyapunov exponent can be calculated by using an effective first order Langevin equation; we also 
recall the exact results for white noise. In section 3, we rederive the rigorous functional evolution equation of P.D.F. ; 
although this equation is purely formal and is not closed (it involves a hierarchy of correlation functions), it will 
be used as a systematic basis for various approximations; we also carry out an exact perturbative expansion of the 
Lyapunov exponent in the small noise limit. In section 4, we consider a mean-field type approximation known as the 
'decoupling Ansatz' which provides a simple expression for the colored noise Lyapunov exponent in terms of the white 
noise Lyapunov exponent. In section 5, we consider two small correlation time approximations that both lead to an 
effective Markovian evolution : we show that these approximations are fairly accurate in the small noise regime. In 
section 6, we investigate the large correlation time limit by performing an adiabatic elimination : this approximation 
is quite suitable for the large noise regime. The last section is devoted to a synthesis and a discussion of our results. 



II. GENERAL RESULTS 



A. The random harmonic oscillator and the Lyapunov exponent 



A harmonic oscillator with a randomly varying frequency can be described by the following equation 

^+7^ + (-^ + eoW)x = 0, (1) 

where x{t) is the position of the oscillator at time t, 7 the (positive) friction coefficient and uj the mean value of 
the frequency. We assume that the frequency fluctuations ^(t) are modelised by an Ornstein-Uhlenbeck process of 
amplitude Vq and of correlation time tq. The long time behavior of x{t) is characterized by the Lyapunov exponent 
defined as 

A(w,7,I?o,ro) - lim ;^(log(^ + ^)) - lim ;^(logi?) , (2) 

t — >oo Zt Z Z t — ^cxj Zt 

where the brackets () indicate an averaging over realizations of the noise between and t, i.e., an averaging with 
respect to the Probability Distribution Function (P.D.F.) Pt{x,x); the quantity E is the energy of the system. 
Taking the time unit to be uj~^, we obtain the following dimensionless parameters, 

a ^ — , Vi ^ — , Ti ^ ujTo . (3) 

UJ UJ-' 

In terms of these parameters, equation becomes 



^ + «^ + (i + e(0)^ = o. (4) 



The Ornstein-Uhlenbeck noise ^(t) now has an amplitude Vi and a correlation time ri and can be generated from 
the following linear stochastic differential equation: 

+ (5) 

dr Ti Ti 

ri{t) being a Gaussian white noise of zero mean value and of amplitude Vi. In the stationary limit, ^(t) has exponen- 
tially decaying time correlations: 

mm) = g^exp(-|t-t'|/ri) . (6) 



When Ti 0, the process ^{t) becomes identical to the white noise. In terms of the dimensionless parameters, the 
Lyapunov exponent is given by 



A([j,7,2?o,T-o) = LjA{a,Vi,Ti) . 



(7) 



The origin (x = x = 0) is a fixed point of equation I^J and in the absence of noise, it is a stable and global attractor. 
If the noise is sufficiently strong, the origin becomes unstable and the system exhibits an oscillatory behavior in the 
stationary state. Thus, the system can undergo a noise induced phase transition from an absorbing state to a non- 
equilibrium steady state (NESS). The Lyapunov exponent vanishes on the transition line between the two phases : 
when A < 0, the origin is stable and when A > 0, the stationary state is extended. In the extended phase, the system 
undergoes random oscillations with increasing amplitude; nonlinearities must therefore be taken into account 23, 24] . 
In the vicinity of the phase transition line, on-off intermittency is displayed and the moments of cc exhibit a multifractal 
scaling behavior |25'|. Thus, for the random phase oscillator, the sign of the Lyapunov exponent determines the phase 
of the system. In figure we plot this transition line computed numerically for r = 1.0 and 2.0. We also draw the 
transition line for the white noise, that was calculated analytically in p^ . 
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FIG. 1: Critical curves obtained by simulating Eq. The curves with black circles (•) and with squares (□) correspond to 
T = 1 and T = 2, respectively. The full black line is the analytical result for the white noise case. 

The objective of this work is to study various analytical approximations for the Lyapunov exponent when the noise 
is an Ornstein-Uhlenbeck process. These approximations will allow us to deduce the analytical features of the phase 
diagram of a random frequency oscillator subject to an Ornstein-Uhlenbeck noise. 

B. Reduced first order equation for the Lyapunov exponent 

As in the white noise case, the Lyapunov exponent can be calculated by solving the following first order nonlinear 
stochastic differential equation 



z = fiz)+m, 
with f{z) = e — , 

and e = sign(a — 2) ; 



(8) 
(9) 
(10) 



the auxiliary variable z is related to the original variable xhy z = ^ + ^. The equation (jSj) is derived from equation Q 
in Appendix 1X1 The noise ^(i) in this equation is an Ornstein-Uhlenbeck process of amplitude A and correlation time 
T, ie., ^(t) satisfies the relation 



1 1 

T T 



(11) 



where ri{t) is a Gaussian white noise of zero mean value and of amphtude A. Thus, in the stationary hmit we have 
{£,it)£,{t')) = ^ exp {—\t — t'l/r) . The parameters e, A and r that appear in the reduced problem are related to the 
dimensionless parameters {a,I)i,Ti) as follows (see Appendix IXjl 



Vi / q;2 
For a < 2 (underdamped case) : e = —1, A = —j- , t = ti \ \ . (12) 

1 /3 

For a = 2 (critical damping) : e = 0, A = l,r = T)^ ti , (13) 
or, equivalently, e = 0, A = ViTi , t = 1 . (14) 

For a > 2 (overdamped case) : e = +1, A = — ^ jj^ ~ '''^ \ ^ ^ ' ^^"^^ 

We show in Appendix A that the calculation of the Lyapunov exponent A(a, ri) of the random oscillator driven 
by noise can be reduced to the calculation of the following quantity : 

A(e,A,r)= mxL (^(z), + i|(log(z2 + l)),) . (16) 

Hereafter, A(e, A, r) will be called the reduced Lyapunov exponent. Equation 116|l seems to involve the time dependent 
P.D.F. of z, but in Appendix B, we prove that 

A(e,A,r) = (z), (17) 

where the stationary average of z is taken in the sense of principal parts. The formula (|17ll requires only the knowledge 
of the stationary P.D.F. of the random variable z satisfying the Langevin equation 0. 

The Lyapunov exponent A{a, ri) of the random oscillator is related to the reduced Lyapunov exponent A(e, A, r) 
as follows (see Appendix A) 



For a < 2, A(a,I?i,Ti) = ^ 1 - — A(e = -1, A, r) - - . (18) 

Fora = 2, A(2, X>i, ti) = X>* A(e = 0, 1, Pj^Vi) - 1 , (19) 

or, equivalently, A(2, X>i, n) = ;^A(e 0, Pirf , 1) - 1 . (20) 



Fora>2, A(a,I?i,Ti) = y - 1 A(e = 1, A, r) - | . (21) 

To summarize, the mathematical problem we have to solve is to find the stationary P.D.F. corresponding to the 
Langevin equation (jHl and calculate from it the first moment of the random variable z. Then, using equations l|7|l. 
(|17|l and H18lll9l or l20ll21|) . we can calculate the Lyapunov exponent in terms of the initial parameters of the stochastic 
oscillator. 

C. The Lyapunov exponent for white noise 

When the noise is white, the Fokker-Planck equation for Pt{z), 

- -Yz - ' + ' 

can be solved exactly in the stationary limit 18] and we obtain 

P(^) = ^exp(^|(ez-y)^y duexp(^-|(eM-y)^ , (23) 

where the current J is determined by the normalization condition dzP{z) — 1. In jl^ . we calculated (z) and 
deduced the following expression for the white noise Lyapunov exponent: 



1 I J+°^duV^e-^(^'-"^^''+'^) 



Awh.te(^^ V,) = -{ :iS> ^ , ... .,3^ -c^} ■ (24) 
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r+oo du „-w((i-'t-)"+t|-) 
Jo 



For small values of the noise amplitude and for a < 2 the Lyapunov exponent admits an expansion in powers of Vi 
that begins as follows: 



A-^^-^a^V,) = ^^^^ - I + OiVl) .^1-^^-1 + O(A^) , (25) 

where A is defined in equation p2|l . The stability boundary of the trivial solution x = i = is given by A^'^'*°(a, = 
0: this equation defines the critical line I'J(a) in the parameter plane (a, Vi). In [l^. we obtained the implicit equation 
for this critical curve and showed that for small and large values of a, we have, respectively, 

and VI ~ Ka^ with K ~ 3.54 . (26) 
At critical damping a = 2, the expression of the Lyapunov exponent is particularly simple 

A^'^"°(2,I?i) = T^('^V-l. (27) 
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III. PROPERTIES OF THE COLORED NOISE P.D.F. 



When the noise £^{t) is an Ornstein-Uhlenbeck process, the stationary P.D.F. of the random variable z driven by 
equation (jHJ cannot be exactly determined. Thus in contrast with the white noise case, there exists no closed formula 
for the Lyapunov exponent A. In the first subsection, we derive, following [2^, the exact functional evolution equation 
for the P.D.F. with colored noise. Although this equation can not be solved, it will be useful in the sequel to construct 
various Fokker-Planck type approximations 26, 27, 28, 29]. In section UlI Bl we derive a perturbative expansion of A 
in terms of the noise amplitude. This expansion will be used in the following sections to test the accuracy of some 
commonly used approximations. 



A. Functional evolution equation for the P.D.F. 



When the noise has non- vanishing time correlations, the dynamics of z is non-Markovian. The evolution of the P.D.F. 
Pt{z) can no more be described by a closed Fokker-Planck type equation. The P.D.F., rather, evolves according to an 
integro-differential equation that involves a memory kernel, non-local in time. Using a functional calculus approach, 
the evolution equation of the P.D.F., defined as, 

Pt{z) ^ {S{z{t) - z)) , (28) 

can be derived with the help of the Furutsu-Novikov formula and is given by 



^ ' - if{^)Pt{z)) + ^9^2 7^ ds exp (-|t - s|/t) ^,5(z(t) - 



dt d 



(29) 



where 5z(t) / 5S^{s) represents the functional derivative of the solution z{t) at time t with respect to the value ^(s) of 
the stochastic process at time s. From equation |(SJ), we find that 



6z{t) 



= 0(i-s)exp^*/'( 



5£,{s) 

where 0{t — s) is the Heaviside function. Substituting this formula in equation (|29|l leads to 
dPt{z) 



(z(w))du (30) 



dt 



d A d^ f* 

-^(/(^mw)+^^y/-xp, 



(-|t-s|/r)^(z(t)-z)exp f'(z(u))du . (31) 



Although this evolution equation for the P.D.F. looks like a Fokker-Planck equation, it is not a closed equation because 
it involves higher order correlations: the functional derivative involves the knowledge of the function z at different 
times. It is only in the case of white noise that equation H31|l reduces to the usual Fokker-Planck equation. Thus, in 
order to make some progress, some closure assumptions must be made. 



B. Small noise perturbative expansion of the P.D.F. 



A closed Fokker-Planck equation associated with the Langevin equation ((SJ can be constructed by a Markovian 
embedding of the coupled stochastic equations iP^andlTT]!. However, we then have to study a two variables Fokker- 
Planck equation for the joint P.D.F. Pt{z,^) : 

This equation does not appear to be solvable even in the stationary limit but it can be used to obtain exact perturbative 
results 30] . We derive here a second order expansion of the Lyapunov exponent for small noise and deduce from it 
the equation of the critical curve near the origin (both Vi and a are small) : we thus take e = — 1. We write the 
perturbative expansion of the stationary solution of equation (|32(l in the form: 



with 



A{z, ^) = 1 + ^ai{z) + ea2{z) + e^zi^z) + . . . , (34) 
B{z,0 = h^{z) + ^bi{z) + . . . . (35) 

The functions and hi satisfy a hierarchy of differential equations of first order in z that can be solved recursively. 
The reduced Lyapunov exponent A(e = —1, A, r) is calculated from the following expression (which has the advantage 
of converging faster than equation H17|l ) 

A(. = -l,A,r) = -^^^, (36) 

where the average is calculated with respect to the stationary measure P{z, ^). Equation H36|) is obtained by multiply- 
ing equation (|32|l by the factor 1/2 log (z^ -I- l), integrating the right hand side by parts and taking the limit i — ^ oo. 
Taking into account the fact that ^ A we observe from equation H36|l that we must determine the functions ai, a2, 
as, 6o and bi if we want to obtain A(e = —1, A. r) up to the second order in A. After some lengthy calculations, we 
find 

, ^ , A A2 r (624r6 1632t5 + 176r4 - 312t3 -I- 137r2 - 324r + 45) ,^,^3, 
^ ^' ^) - WT^) + 32(1 + .^)(l + 4.^)3(9 + 4.^) + ) • 

This expansion provides the local behavior of the critical curve near the origin. Using equations H12l and [T^ and 
retaining only the terms of the first order in a, we find the following equation for the critical curve 'D1{o) : 

4(1 + 4.^). . VI + (^^)- r + 1632r^ + 176r^ - 312.3 ^ ,3,,. _ 3,4, ^ 43) 

^ ^ ^ 4(l + r2)(l +4T2)2(9 + 4r2) ^ ' 

IV. THE DECOUPLING ANSATZ 

The decoupling Ansatz can be seen as a mean-field approach in which higher order correlations are neglected. 
As such, it is the simplest approximation from the technical point of view. It has the advantage that it does not 
require any a priori assumption on the correlation time of the noise. Starting from equation (|29|l . we assume that the 
expectation value factorizes as follows : 

S{z{t) - z)^'j ^ (5{zit) - z)) ^ e{t - s)P,{x) ^exp ^* /'(z(«))du) ^ , (39) 

where the last equality is obtained using equations H28land l5mi . Using again a decoupling approximation and taking 
the stationary limit t,s ^ 00, wc find 

^exp ^* /'(z(u))du^ ^ = exp f'{z{u))d^^ = exp {{t - s){f)) , (40) 



where the mean- value (/') is calculated in the stationary state. After inserting these approximations in equation H31() 
and calculating explicitly the integral obtained, we find the effective Fokker-Planck equation for the decoupling 
approximation: 



A 



(41) 



When T ^ this equation becomes identical to the white noise Fokker-Planck equation. As usual in mean-field 
approximations, the solution of equation (|41|) is obtained by solving the non-correlated case (i.e., the white noise 
problem) supplemented with a self-consistent condition. Using the expression ^ for the function /, we obtain the 
following decoupled effective Fokker-Planck equation 



dPt{z) _ d , A d^P,{z) 

~dr - 'dl ' ' >^'^'^> + 2(1 + 2r(^)) dz^ 



(42) 



This equation is identical to equation H22(l obtained for white noise, except for the noise amplitude A that is replaced 
here by A/(l + 2t{z)). Recalling from equation p7|l that A(e, A,t) — (z), we obtain the following relation 



A(e,A,r)=A 



white 



A 



_H-2rA(A,r) 

Reverting to the original parameters with the help of equations (|18I - I21|) . we conclude that 



A(a,I?i,ri) = A 



white 



a. 



l + ari + 2riA(a,Di,ri) 



(43) 



(44) 



with A*'^'*'' given by equation (|24|l . Equation H44|l is an implicit expression for the Lyapunov exponent of the harmonic 
oscillator subject to Ornstein-Uhlenbeck noise (this is an approximation because it is derived from the decoupling 
Ansatz). The critical line I?i(a) is obtained by taking A(a,I'i,Ti) = and its equation is 



A(a,I?i,ri) = A 



white 



a. 



1 -I- aTi 



(45) 



The critical line for the Ornstein-Uhlenbeck noise, as obtained from the decoupling Ansatz, is thus readily deduced 
from the white noise critical line by a simple coordinate transformation. The resulting phase diagram is drawn in 
figure 121 it agrees qualitatively with the numerically computed curves given in figure^ 
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FIG. 2: Phase diagram obtained by the decoupling Ansatz for r = 0.1, 0.5, 1.0 and 2.0. The full black line corresponds to the 
white noise problem. 



Using the expansions (|26() obtained for the white noise, we derive the asymptotic behavior of the critical hne in the 
decouphng approximation : 



VI ~ Aa for a < 1 (46) 
VI ~ 3.54Ta^ for a > 1 . (47) 

For small values of the dissipation rate a, the critical curve obtained by the decoupling approximation is linear with 
the same slope as for white noise : comparing equation H46|l with the exact expansion H38|l , we observe that this result 
is quantitatively correct only if r ^ 1. For large values of the dissipation rate a, equation (|47|l shows that Vi scales 
as the fourth power of a; we recall that for white noise, the asymptotic scaling is V^ cx (see equation (|26|l l. This 
asymptotic scaling, different from the one obtained in the white noise case, is valid as soon as r 0, however small 
T may be. This feature will be investigated more specifically in section 6. 

Finally, we can study the special case a ~ 2. Defining L = ti ^A(2, ti) + 1^7 we find from equations H27|) and 

(03 that L satisfies 

3 7r3/2 

^ ^^ + ^^^-TfxW ^^^^ 

Solving this quartic equation for L leads to the Lyapunov exponent for the critical value a = 2. In particular, we have 

M2.V,.r,) . (^)*(|^)' -1.0.331 (^)'-l fa I>,.f»l. (50, 

The polynomial equation (|48(l interpolates between the white noise solution (equation 1)49(1 is identical to equation ((27|l l 
and the colored noise scaling with exponent 1/4 when Pirf ^ 1 (equation 1)50(1 ). 



V. EFFECTIVE FOKKER-PLANCK EQUATIONS FOR SMALL CORRELATION TIME 

In this section, we study the effective Fokker-Planck equations that are valid when the correlation time r of the noise 
is small. The stationary solutions of these_ equations allow us to calculate small correlation time approximations of 
the Lyapunov exponent. In the spirit of ^2^, we derive these approximations from the exact evolution equation lp?T|l . 



A. First order effective Folcker-Planclc equation 

In the small correlation time approximation, we assume that r <C 1 and thus, in equation 1(31(1 . we have |i — s| 1- 
We then make the following approximation: 

exp ^* nz{u))dv}j c exp {(t - s)f'{z{t))) . (51) 

Inserting this expression in the functional equation 1(31(1 and using equation ((28(1 . we obtain an effective Fokker-Planck 
equation, valid at first order in r 



dt 



- {J{z)Pt{z)) + - — ((1 + r/'(z))P,(z)) . (52) 



We remark that in the limit r ^ this equation reduces to the white noise Fokker-Planck equation. Using the 
expression (0 for the function /, we find that the first order effective Fokker-Planck equation leading to the Lyapunov 
exponent is given by 

2^ = _|(,._,V,,.,)l-|i^(,l-2„,P.w). ,53, 



We now solve this effective Fokker-Planck equation in the stationary hmit. Introducing the stationary current J we 
obtain 



J = (z^ - e)P{z) + ^1: ( (1 - 2rz)P(z) ) . (54) 



2 dz 

This equation is solved by the variations of constants method. In terms of the parameters 

20 ^_ and ^=_^, (55) 

the stationary P.D.F. can be written as 

2J, , / z2 z \ exp(-2^-2^ 

Foi z<zo, Piz) ^ —\2tz -l\^-'exp (— + ——] du ^ 



A' ' '^\2tA 2t^A / _ |2rw-l 



2J, 



For . > zo, Piz) . ~-i2rz - 1)^- exp [^ + ^)J^ ^"^27^1^1)^ ' (5^) 

We remark that when r ^ this P.D.F. becomes identical to the formula 123() derived for white noise. Using 
equation we can show that P{z) ~ -p- when z ^ ±oo . The behavior of P{z) when z ^ zq is as follows: 

For A<1, P{z)^- (57) 

For P(z) - ln|z - zol , (58) 

For A>1, P{z,) = -^±-^^. (59) 

Thus, for A > 0, the stationary P.D.F. is a positive and normalizable function and the current J is fixed by imposing 
J-^ P ~ 1- The small correlation time approximation breaks down when A < 0, i.e., when 4T^e > 1; this happens 
in the case e = +1 and r > 1/2. In the rest of this section, we consider only the case r < 1/2 : the approximation is 
then well defined and can be used to study the Lyapunov exponent. 

For small noise, A ^ 1 and low damping, e = 1, we derive a perturbative expansion of the stationary P.D.F. (|5(: 



A(. = -1, A, r) = I + ^ + 0(A3) . (60) 

We observe that the linear term in A agrees with the exact result (|37|l if terms of order ©(r^) are neglected. However, 
the term proportional to A^ does not coincide with the exact result even at the first order in r. 

For large noise, A ^ 1, and for a given correlation time r, equation (|55|) implies that A ^ 0. From equation H57|l . 
we conclude that the stationary P.D.F. becomes more and more localized in the vicinity of zq as A grows : hence 
(z) zq = 1/2t as a ^ oo. Thus, the first order effective Fokker-Planck equation predicts that the Lyapunov 
exponent saturates to a finite value when the amplitude of the noise grows. This prediction is unphysical and in 
contradiction with numerical results. The first order effective Fokker-Planck equation is therefore meaningful only for 
small amplitudes of the noise. 

In conclusion, the first order effective Fokker-Planck equation does not provide a good approximation to calculate 
the Lyapunov exponent. In the next section, we discuss a more sophisticated approach that leads to good results for 
the small noise regime at least. 



B. The 'Best Fokker-Planck Equation' 

This approximation has been first proposed in |3lj |. Starting from an approximate integro-differential equation 
for the P.D.F. valid at first order in the noise amplitude, the evolution kernel is calculated by a resummation to all 
orders in the noise correlation time. This procedure results in an effective Fokker-Planck equation, called the Best 
Fokker-Planck Equation (B.F.P.E.). Although this equation is not free from inconsistencies [S^], it provides in some 
cases useful insights that agree qualitatively with numerical simulations. (Another approach, that we shall not follow 
here, is to reject the assumption of stationarity and to study non-conventional diffusion regimes of a Fokker-Planck 
equation with a time-dependent diffusion constant H^.) 



Following pq . we derive the B.F.P.E. approximation from the exact functional evolution equation H31() . In the 
B.F.P.E. approach the amplitude A of the noise is supposed to be small and contributions of order A" with n > 2 
are neglected. In order words, the exact evolution equation H31|l is replaced by 



dPtiz) 



{f{z)Ptiz)) + 1^ ds exp (-|t - s|/r) exp (^J^ f'(z(u))du 



(J(z(t)-z)), (61) 



dt 

where z{u) is the solution of the deterministic ( i.e., noiseless) equation 

dz 

— = f(z{u)) with z(<) = z. (62) 
du 

Substituting this equation for z in equation H61|l . we derive the B.F.P.E. 

^ -I 4l^(-'<--)«'»)) "'^'> I! '-''^'^m^^ 

The space and time dependent diffusion factor Dt{z) can be explicitly evaluated for f{z) — e — z^ . In the long time 
limit t — > 00, we obtain 34] 



dt dz 



if{z)Ptiz)) + ^^^{Diz)P,iz)) with i?(.) = . (64) 



The effective diffusion coefficient D{z) is everywhere positive when e = —1 or for all values of t. When e — 1, D{z) 
is everywhere positive only if r < 1/2. In all these cases, we obtain 

P{z) = — exp (-4tAz) D{z)-'^^-'^ / duexp (4tAu) D{uf^ , (65) 

where A was defined in equation 1)55(1 . The value of J is again fixed by normalizing P. Again, when t ^ 0, this 
expression becomes identical to the formula H23|l for the P.D.F. obtained in the white noise case. From this expression 
of the stationary P.D.F. we derive the expression of the Lyapunov exponent: 

_ JZ,^ dz z exp i^ArAz) D{z)-^^-' Au exp (4rAu) D{uY^ 
'^^ /^^dzexp(-4TAz) i?(z)-2^-i/f^duexp(4rAw)i^(u)2^ ' ^ ^ 

where the integral in the numerator is defined in the sense of principal parts. This closed formula can be used for a 
numerical evaluation of the Lyapunov exponent in the B.F.P.E. approximation (see figureEJ. 
The perturbative expansion of the Lyapunov exponent (|65|l when A <C 1 and e = 1 is given by 

Comparing this expression with the exact result ((37|l . we observe that the first term of this expansion with respect 
to A is correct. The B.F.P.E. yields the exact result to all orders in r and performs indeed a complete resummation 
with respect to the correlation time. However, the term in A^ does not agree with the exact result, even in the 
r ^ limit. We recall that the BFPE approximation is intrinsically a small noise expansion and we see clearly, in 
this specific example, that this approximation is not valid at second order in the noise amplitude. This agrees with 
the general analysis of that the BFPE can not be used for moderate values of the noise because the neglected 
non-Fokker Planck terms have a contribution of the same order as the terms obtained after resummation with respect 
to the correlation time. 



VI. ADIABATIC LIMIT FOR LARGE CORRELATION TIME 



In the previous section, we studied small correlation time effective equations which yield fairly good approximations 
for the critical line for small-to-moderate values of the noise amplitude. But these approximations always break down 
when a ~ 1/r: a specific approximation scheme is therefore needed for this range and will become more and more 
relevant as r grows. 



In this last section, we thus study the large r case. When the correlation time is large, the Ornstein-Uhlenbeck noise 
becomes a slow variable and the dynamics of the random variable z can be simplified. We shall first perform an ele- 
mentary adiabatic elimination that yields an expression of the Lyapunov exponent in the r — > oo limit. Then, we show 
that these results can be obtained in a more systematic manner thanks to the Unified Colored Noise Approximation 
(UCNA). 

When r is large, the noise ^(t) is slowly varying and we make the approximation that z keeps adjusting itself to 
the stationary solution of equation 0. We consider here only the domain a > 2 (i.e., e = 1) for which the small 
correlation time expansions do not provide reliable results. Thus, we study the limiting case where the noise f is 
quenched; the equation ((SJ then reduces to 

i = (68) 

where ^ is a time-independent Gaussian random variable of variance A/(2r). The value of z is then given by the 
stable fixed point of equation Ht)8|l if such a solution exists. We must distinguish two cases: 



• if ^ > —1, then z = + 1 (the solution z = — V^~i-T is unstable); 

• if ^ < —1, then equation H68() has no fixed point. The random variable z is given by z{t) — ylcotan(Ai) with 

= —{^ + 1). This running solution implies that z is distributed over the whole real axis. 

From this discussion, we conclude that the distribution of z is given by 



For z > 0, P{z) = 2^/ — z exp (-r(z^ - if /A) + J -r I ^ , exp (-rf /A) dC , 




for z < 0, P(z) ^\—r \ \, exp {-re /A) d^ . (69) 

V 7I-A 7r(z^ - C - 1) 

This P.D.F. consists of two parts: one term is of Lorentzian type and is an even function of z; the other term exists 
only for z > and is a rapidly decaying function. This contribution provides a strictly positive value for the mean 
value of z. Using equation (|69|l . we calculate the Lyapunov exponent 

A(e = l,A,r) = (z) =2,/^ /"z2exp(_^(^2_^)2/^)j^^ IW4) /Ay/4 ^^^^^ A ^ ^ ^^^^ 
V ttA Jq 2v7r V r / r 

This simple adiabatic approximation provides the behavior of the Lyapunov exponent for large amplitude of the noise 
and it predicts a scaling in agreement with the decoupling Ansatz (see equation ()47|l'). Moreover, it explains the 
existence of long tails for the P.D.F. of z near infinity (P(z) ~ 1/z^ when |z| oo) and the asymmetry in P(z) for 
z —z. Because of this asymmetry that favors positive values of z, the Lyapunov exponent, (z), is always strictly 
positive. Using equation H7()|l. we obtain the asymptotic behavior of the critical curve for large values of the noise 
amplitude, 

This approximation is fairly accurate even from quantitative point of view and provides a good approximation of the 
critical noise amplitude even at moderate values: this can be seen in figure |21 where the expression H71(l is compared 
with the critical curve, computed numerically. Moreover, the asymptotic behavior in equation H71(l agrees well with 
the prediction of the decoupling Ansatz, equation H47|l . 

For the special case a = 2, we can repeat the above calculations and find that 

This expression is very similar to equation H5()|l obtained by decoupling Ansatz. However the limit t — > is ill-defined 
here; this is not a surprise because the adiabatic elimination makes sense only if r is large. 

Finally, we show that the adiabatic elimination can be performed in a systematic manner by using the simplest 
version of the unified colored noise approximation |2^ 123 . Is^ l (for more sophisticated UCNA schemes, see e.g., 



[3a,|33). Taking the time derivative of equation © and using equation (fTT|l . we obtain 



/'(z)z + e-/'(^)i--e+-ry. (73) 

T T 
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FIG. 3: Various approximations for the phase digram at r = 1. The critical curve, computed numerically, is plotted in full line 
with black circles (•). The small and large r expansions, equations 1)37^ and H71|l respectively, are drawn with dashed lines. 
The decoupling approximation, equation 1451 , is plotted with a dot-dashed line. The BFPE approximation is represented with 
A. 



Using again equation (jSl to express the Ornstein-Uhlenbeck noise in terms of z and the white noise ry, we derive the 
following second-order nonhnear Langevin equation driven by white noise: 

(Pz dz 1 1 

^+7rW^-/(^) = ^r;(T) where ^,{z) = -j= - ^f'iz) , (74) 

where the new scale of time is T = t/^. The effective damping coefRcient 7r(-z) diverges when r ^ and r — > oo. 
After an adiabatic chmination of the incrtial term (vahd for 2tz + 1 ^ ^/t) the foUowing effective stochastic equation 
is obtained: 

TF^^—r-{+ lu ( M T), with j,{z) = ^+2V^z. (75) 

The solution of the associated Fokker-Planck equation leads to the stationary P.D.F., which in the domain of validity 
of the UCNA is given by 

P(z)=7V|2Tz + l|exp-— — + — -erz^-ez for z > . (76) 

A y 2 o J ZT 

The constant M is fixed by the normalization condition on P{z). From this expression we deduce that the Lyapunov 
exponent is given by 

z\2tz + 1| exp-^ ("^ + T ^ '^'''^'^ ^ '^A 

A(e=l,A,r) = — — — ^ }rz^ , .3 ^ (— ■ (77) 

Jp |2tz -I- 1| exp-^ (^ + ^ - eTz2 - dz 

Here, the range of the variable z is taken from to cxd instead of z > — ^ (this is a good approximation when r is 
large). For e — \ and A 3> 1, this expression scales as (A/r)^/'', in agreement with the behavior predicted by the 
decoupling Ansatz and the adiabatic elimination. For the special case a = 2, equations l|2()l and[77 | lead to 

A(2,I?i,ri) = , - l ^ -^,7:^7:^^ l ^ l - 1 for V,tI « 1 , (78) 



1 /o°°^|2^ + l|exp-p^(^ + tjdz 1 fWi 

n /;°|2z-|-l|exp-^(^ + f)dz ' " r(2/3) 



2 



(79) 



Thus, the UCNA approximation permits an interpolation between small and large values of the correlation time. 
The predicted behavior matches the white noise scaling for small t, and the colored noise scaling at large r. We remark 
that prefactors in equations 1)781 and [7(^ are different from those obtained in equations H27|l and (|7()|l . respectively; 
this is due to the fact that in the UCNA the stationary current vanishes : the running solution of z are overlooked 
and therefore this approximation does not describe satisfactorily the tails of the P.D.F. of z. 

This study of the r ^ oo limit shows that the absorbing phase becomes more and more stable as the correlation 
time increases. Besides, it confirms that the existence of a time correlation of the noise renders the critical curve 
steeper and modifies its scaling behavior for large values of the noise amplitude. 



The long time behavior of the stochastic oscillator with random frequency is controlled by the sign of the Lyapunov 
exponent. In the case of a white noise perturbation, this exponent can be calculated exactly and the phase diagram 
of the stochastic oscillator can then be rigorously determined. The aim of this work is to study the effects of time 
correlations on the noise-induced bifurcation of the stochastic oscillator. In the case of an Ornstein-Uhlenbeck noise, 
an exact calculation seems to be out of reach and therefore we use various approximation schemes in order to derive 
analytical expressions for the Lyapunov exponent and to draw the phase diagram. Since the different approximations 
have distinct regions of validity, their study has allowed us to derive a global picture of the behavior of the system 
in the parameter space. In particular, we have derived the scaling behavior of the phase boundary in regions where 
the amplitude of the noise is small or large. Our results agree fairly well with numerical simulations and with exact 
perturbative expansions. These comparisons allow us to test the validity of the different approximation schemes. 

We remark that the effective first order equation (jSJ, that we have used to calculate the Lyapunov exponent, can be 
interpreted as describing an overdamped Brownian particle driven by colored noise in a metastable potential. For such 
a system, a fundamental and extensively studied quantity is the mean escape-time of the particle jlj from the well. 
The relation between the Lyapunov exponent and this mean escape-time deserves to be clarified. In particular, one 
could then use path-integral techniques JsS, to study the phase diagram of the random frequency oscillator. The 
variation of the Lyapunov exponent with the correlation time of the noise could then be related to the phenomenon 
of noise enhanced stability in metastable states |40| . 

The study carried out here for the oscillator with multiplicative noise can be adapted to other stochastic systems. 
For example, Schimansky-Geier et al. |4ll | have shown that the random Duffing oscillator with additive noise undergoes 
a subtle phase transition that does not manifest itself in the stationary P.D.F. (which is simply given by the Gibbs- 
Boltzmann formula) but affects the properties of the random attractor in phase space. This problem is mathematically 
equivalent to the noise-induced bifurcation of a linear oscillator subject to a multiplicative colored noise. This noise 
has a finite correlation time and if we approximate it by an Ornstein-Uhlenbeck process, the system becomes identical 
to the one studied here. 



It has been suggested recently [42j that the Poisson process may be suitable as a paradigm study of time correlation 
effects in random dynamical systems. In fact, in many cases, a complete analytical study can be carried out for a 
random variable driven by Poisson noise, but not when the stochastic force is an Ornstein-Uhlenbeck process. We 
have indeed carried out the exact calculation of the Lyapunov exponent of a random frequency oscillator subject to a 
Poisson noise and obtained the exact phase diagram of the system ji^l • The results are in qualitative agreement with 
those found here with an Ornstein-Uhlenbeck noise. 



In this appendix, we show that the calculation of the Lyapunov exponent can be reduced to solving a nonlinear 
ffi'st order Langevin equation. We introduce the function y{t) — ^, that satisfies 



VII. CONCLUSION 



APPENDIX A: DERIVATION OF THE EFFECTIVE LANGEVIN EQUATION 



dy 



(Al) 



dt 



Using the fact that 




(A2) 



■ 2 2 

where £^ = ^ + ^ is the total energy of the system, we derive the following identity: 




(A3) 



This identity implies that the Lyapunov exponent is given by 



A(a,I?i,Ti) = i hm ^(logii;) = ^Um (^{y)t + ~ {\og{y^ + 1))^^ . (A4) 



2 t^oo dt 

We introduce an auxihary variable z defined as 



Eliminating y from equation (|A1(I . we obtain 



X a a , 



z^^-l-z^ + m- (A6) 



We now show that the dissipation parameter can be eliminated from this equation by a suitable redefinition of the 
parameters involved in the problem. We have to distinguish three cases: 



• Underdamped case (a < 2): in terms of i y 1 — ^ t, the evolution equation of z becomes 

z^-l-z^+m. (A7) 
where ^ is an Ornstein-Uhlenbeck noise of amplitude and correlation time given by 



• Critical damping [a — 2): the evolution of z is given by 

i = + (A9) 

1 /3 

We can rescale the time variable either as t :— t, or as t := t/ri. The amplitude and correlation time of ^ 
are then given, respectively, by 

A =1, T = vl^^Ti. (AlO) 

or A = Pirf , T = 1 . (All) 

We notice that in the critical damping case there is only one free parameter in the problem (in the white noise 
limit Ti = 0, no free parameter is left). 



Overdamped case (a > 2): in terms of t :— y ^ — I t, the evolution equation of z becomes 

z^l- z^ +£.{t), (A12) 
where the amplitude and correlation time of ^ are given by 



We have thus shown that the calculation of the Lyapunov exponent of the linear oscillator with random frequency 
can be reduced to the study of the following equation 

i = e-z2 + ^(<), (A14) 

with 



e = sign(a - 2) , 



(A15) 



i.e., e = —1,0 or 1 when a < 2, a — 2 and a > 2 respectively. Using equation (|A4|) and taking into account the 
rescahng of time, we conclude that 



A(a,I?i,Ti) = KA(e,A,r)-|, (A16) 



where we have defined 



1 d 



t^oo V ^ 2dt\ V 2 



a 



A(e,A,r)=lim (^), + -- log + l , (A17) 



/ 2 1 /3 / 2 

and the coefhcient k is equal to wl — ^ , or J ^ — 1 when a < 2, a = 2 or a > 2, respectively. We remark 
that when t ^ oo, 

(z-f)^+l 
2' ' 'I/t ' '"^ \"'° z2 + l 



log [{z - '^f + l) )^ - (log + log ^ > , (A18) 



where the average on the right hand side is taken with respect to the stationary P.D.F. of z and does not depend on 
time. Therefore, equation (jA17p is equivalent to 



A(e,A,r).mrr ((.). + l|(log(z^ + l)),) 

which is identical to equation 



(A19) 



APPENDIX B: PROOF OF EQUATION ljrf|l 

We now derive the identity p7|l for the Lyapunov exponent that involves only the stationary P.D.F. of z. In all the 
cases considered in this work, we obtain an effective Fokkcr-Planck equation of the type : 

dPtjz) _ , A & 

dt dz 

In the stationary limit, we have 



(/(z)P,(z)) + |^(i?(z)P,(z)). (Bl) 



J = -f{z)P{z) + ^±[d{z)P{z)) , (B2) 
where J represents the stationary current. After an integration by parts, we deduce from equation (jBip that 

d H , s\ 1'^°° , 2z / , , , A d{D(z)Pt(z))\ , , 

Thus, we have 

(^^^^\ lT^^^ ~fiz)Pt{z) + ^£mz)Pt{z)) ^_^ 

Taking the stationary limit and using equation l|B2p . we find 

Im^ {izh + ^1 (log(.^ + 1)).) ^ . (P(.) - ^) . (B5) 



Combining equations (|A19p and ljB5p . we obtain the following identity for the Lyapunov exponent 



AteA.,)^/;;.(pw-^)d=./; 



zP(z)dz , (B6) 



where the last equality has to be understood in the sense of calculating the 'principal part' of the integral, i.e., 

f+oo p+M 



zP{z)dz^ lim / zP{z)dz. (B7) 



Thanks to the formula ljB6|l . the Lyapunov exponent is expressed in terms of the stationary P.D.F. of z and equa- 
tion l|B6p can be rewritten as 

A=(z), (B8) 

which is identical to equation (|17|) . We emphasize that the first moment of the stationary P.D.F. needs to be defined 
only in the sense of the principal values. 
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